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Abstract 



We use quantum field theory methods to study the network model for the 
spin quantum Hall transition. When the couplings are fine tuned in a certain 
way, the spin and charge degrees of freedom, corresponding to the supercur- 
rent algebras su(2)q and osp(2|2)_2 respectively, decouple in the renormaliza- 
tion group flow. The infrared fixed point of this simpler theory is the coset 
osp(4|4)i / sn(2)o which is closely related to the current algebra osp(2|2)_2 but 
not identical. Some critical exponents are computed and shown to agree with 
the recent predictions based on percolation. 
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I. INTRODUCTION 



Recently a number of works have appeared investigating a new universality class of de- 
localization transition referred to as the spin quantum Hall effect. This transition can occur 
in certain dirty superconductors with unbroken su(2) spin-rotation symmetry. Kagalovsky 
et. al. constructed an su(2) invariant network model for the transition and numerically 
determined some critical exponents Senthil et. al. on the other hand modeled the 
phenomena with a supersymmetric spin chain 0J3[]. Remarkably, the network model was 
mapped directly on the lattice onto classical percolation by Gruzberg et. al. and exact ex- 
ponents were computed which agreed very well with the numerical simulations of the super 
spin chain f|. 

In the quantum field theory approach to derealization transitions, the computation of 
critical exponents is normally a difficult strong-coupling problem requiring the existence of 
non-trivial infrared fixed points. The claim that for the spin quantum Hall effect the infrared 
fixed point of the disorder averaged effective field theory is simply percolation is rather 
unexpected and for this reason we set out to understand this using quantum field theory 
methods. Our starting point is the hamiltonian formulation of the network model given in 
0. We carry out the disorder averaging using the supersymmetric method in conjunction 
with conformal field theory methods. This leads to an effective action which consists of 
a conformal field theory with an osp(4|4) super-current algebra symmetry perturbed by 
certain marginal operators which are bilinear in the supercurrents. 

The one-loop renormalization group (RG) /5eta functions we compute for the three inde- 
pendent couplings appear to be as complicated as those for the usual quantum Hall transition 
|| . There exists a fine-tuning of the network model couplings wherein two of the couplings 
are essentially identified and the resulting model is remarkably simpler and can be solved. 
This is due in part to the fact that the stress tensor for the osp(A\A) current algebra can be 
written as the sum of the stress tensors for the osp(2\2) level —2 and su(2) level current 
algebras corresponding to the charge and spin degrees of freedom respectively. This "spin- 
charge separation" is present in the disorder-averaged effective theory and this leads to two 
decoupled /3eta functions. This simplification allows us to determine the non-trivial infrared 
fixed point: the su{2) level degrees of freedom decouple in the flow and the infrared con- 
formal field theory is simply the coset osp(4|4)i/su(2)o- In this way we recover some of the 
percolation exponents predicted in j4|. 

Due to the logarithmic nature of the above conformal field theories, we find that in 
spite of the separation of the stress tensor into commuting pieces, T osp ( 4 | 4 ) = T osp (2\2) + 
T su {2)i the Hilbert space does not factorize. This leads to the peculiar result that the coset 
osp(4\4:)i/su(2)q is not equivalent to the osp(2|2)_ 2 current algebra theory, even though this 
coset possesses the current algebra as a symmetry and the conformal dimensions of the coset 
are the same as for the osp(2|2)_2 current algebra. 

In the last section we study the possible universality classes based on the 1-loop RG 
equations and suggest that the network model is universally attracted to a so-called "strange 
direction" . 
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II. THE MODELS 



Kagalovsky et. al. gave a hamiltonian formulation of their network model |T|. The result 
is the 4x4 matrix hamiltonian: 

H = (t xPx + T zPy ) ® 1 + 1 <g> a ■ a (2.1) 

where f, a are two copies of the Pauli matrices, p x>y = —id x , y , and a is a random spin poten- 
tial. This hamiltonian has the defining properties to belong the class C of the classification 
introduced in ref. ||. Let us perform a unitary transformation H — > (W ® l)H(U <8> 1), 
where U corresponds to a rotation about the x axis for the Pauli matrices: U*t z U = r y , 
U'TyU = —t z , Wr x U = t x . After including su(2) gauge potentials top, one has the following 
2x2 block structure: 

H=( *'* A ~ id ! + M (2.2) 
\-id z + A z a-a ) y > 

where d z = d x —id y , d- = d x + id y , A z = A x — iA y , A z = A x + iA y . In the above hamiltonian, 
a(x,y) is a real, random spin potential, and A^ = ^ a A a {x,y)cr a are random su{2) gauge 
potentials, with A a x real. In the sequel we will take them to have the following gaussian 
distributions 

d 2 x _ . A 

— — a(x) ■ a(x) (2.3) 
2n J 

d 2 x \ 

— a;(x)a%x)) 

As we will show in the next section, renormalization of the effective action obtained upon 
disorder averaging leads to an additional interaction which can be viewed as arising from a 
random mass m(x,y). The complete hamiltonian which leads to a renormalizable effective 
action is then 

H=( a :? + ™ Z id l +A *) (2.4) 
V — %o z + A z a ■ a — m J 

We will take m(x) to have the gaussian distribution 

P(m) = exp (-— [ — mix) 2 ] (2.5) 

V 9m J 27T J 

Note that since d\ = —d^, A\ = A-, the hamiltonian is hermitian if m, a, A x>y are real. In this 
situation the couplings g a ,9m,9A are the variances of normalizable gaussian distributions if 
they are real and positive. Negative couplings g a ,m,A can be interpreted as corresponding to 
imaginary random potentials; hermitian hamiltonians can then be constructed by doubling 




the number of degrees of freedom, as in [ 10] . 



The single particle Green functions are defined by the functional integral 
Z^ 1 J D^*D^ exp (-S) with Z the partition function and 

S= f^- $*{x)i (H - £) (2.6) 
J 2n 
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where £ = E + is. For e = + , this defines the retarded Green function 

G R (x,x';E) = * Is 1 ) = (2-7) 



*= (1 + ), ** = (V-,tU (2-8) 



Letting 

where ^+ is a 2-component fermion if)\, i = 1, 2, and similarly for ^? ± , one finds 

S= / y^(V'-(^ + iA^ + +^_(a z + iA^ + + ia. (V>-dty+ +^_oty+) 

+ im(^-V' + - -iS$ E ) (2.9) 

where 

(We have suppressed the su(2) indices, i.e. ifj-i> + = HiV'-VV' etc) 

In the perturbation H e — er y ®l was added and critical exponents for e were measured 
numerically. Performing the unitary transformation U defined above, this results in the 
following perturbation of the action: 

S e = ~%e J $e = VM^+ ~ (2-10) 

Non-zero e corresponds to a non-zero average random mass m. 

Let us attempt to compare this with the model considered by Senthil et. al. ||. There, 
one had a one-dimensional lattice in the x-direction with sites labeled by j, and a contin- 
uous y direction. At each site there are fermionic degrees of freedom Xj{y)i where j even 
corresponds to left-movers and j odd to right-movers. The hamiltonian is 

H = S / d v({-) 3 x){y){-id y + Vj(y) ■ o)xj{y) - tt]{y){x)+iXj - xjxj+i) 

3 

+ tj(y) ■ (x]+i°Xj + x)pxj+ii) (2.11) 

Taking a continuum limit, Xj(y) ~^ 4>+( x ,y) f° r 3 even and Xjiy) ~^ 4>+( x ,y) f° r 3 odd, 
Vj(y) -> Vy( x ,y), and tj-(y) -> t(x,y) , one finds 



iif = / dxdy(ip_(d y + ■ a)i/)+ - ip„{d y + iff y ■ a)ip + + t (V'-V'h- _ V'-'0+) 

+ zf • (^-#, + ?_<J^+)) (2.12) 



Rotational invariance of the kinetic terms can be restored by adding d x + iff x ■ a to the 
derivatives. Performing a rotation to euclidean space y — > — «?/, one then obtains the action 
Q2.9| ) with a = £, v4 2 = (ff x — iff y ) ■ a and t = im. 

In 0] , t and to were taken as real gaussian distributed but with the same variance. For 
the model defined by (|2.9| ), this corresponds to imaginary m. Letting m — > —im, one sees 
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that for the effective theory obtained after disorder averaging, imaginary m corresponds to 
negative g m . Since the variances are the same, this corresponds to g a = — g m .[] In summary, 
though we have not described an exact mapping between the network model in |1[ and the 
super spin chain in ||f3|], it appears that they should be related on the line g a + g m = 0. In 
fact, one obtains directly the Heisenberg type superspin chain upon disorder averaging only 
along this line ||. We will provide further support of this statement based on symmetry 
in the sequel. As we will show, this line has some rather special properties which allow the 
model to be solved. 

III. EFFECTIVE ACTION AND RENORMALIZATION GROUP ANALYSIS 

A. Supersymmetric Disorder Averaging 

Since we are dealing with a free field theory, the supersymmetric method for disorder 
averaging can be used. Conformal field theory techniques in conjunction with this method 
was used for other models in ^] 0. We augment the theory with bosonic ghosts (3±, /3±, 
i = 1,2, so that the inverse of the fermionic partition function Z(a,m, A) = J Dipe~ s ^ is 
represented as a bosonic functional integral: 

Z(a, m, A)' 1 = J D(3 e" 5 ^) (3.1) 

One can then perform the gaussian integrals over the random potentials. The result is the 
effective action: 

(Px 

Seff = S cft + [ (9a O a + g m O m + g A A ) (3.2) 
J Ztt 

The conformal field theory S c f t has Virasoro central charge c = and the action: 

3 c ft = / ^ (V-<3^+ + $-0$+ + (3-dz(3 + + P_d z p + ) (3.3) 

The first-order action for the bosonic ghosts of conformal dimension ±1/2 can be treated as 
in ||. The operators which perturb away from the conformal field theory are: 

O m = \ - ^_^ + + P.P + - P_{3 + ) 2 (3.4) 



1 This appears to be inconsistent with hermiticity, since the hamiltonian ( p.4| ) is not hermitian for 
imaginary m, suggesting the two models cannot be naively compared in this way. 
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B. /3eta Functions 



The one-loop renormalization group (3eta functions can be deduced from the operator 
product expansions of the marginal perturbing operators Cj ||. Namely, if 

0^)0,(0) ~ O fc (0) + .... (3.5) 

zz J 

then to lowest order the /3eta function /3k = dgk/dlogl, where I is a length scale, is given by 

A = -E c S^i ( 3 - 6 ) 

We normalize the Pauli matrices as follows: o ■ a = (cr 3 ) 2 + o~ + a~~ + o~~cr + , with 
1 \ . /-/0 \\ _ /r/0 



^ • -~ = V2 (3.7) 



<T> " 1 o -\) ' a v " vo oy ' v " vi o 

With this normalization, one has Tra a a b = 25 ab , and f abc f abd = —S5 cd , where [a a , o b ] = 
f abc a c . One also needs the operator products: 

V+(*ty-(0) ~ V-(*)^+(0) ~ 1/2 (3.8) 
^+(2)^(0) ~ -^(^+(0)^1/2 

and similarly for the right-movers ■?/>, /3 with z replaced by z. One then finds 

O A (2,z)O A (0) ~ ^O a , O m (z,z)O m (0) ~ 4^ (3.9) 

ZZ 

a (2,l)O a (O)~4(2Oa-8O A ) 

zz 

1 / 1 „ 3 



A (z,z)0 Q (O) ~ — (--0 Q - -O m ) 

o A (z,2)o m (o) ~ ^ (-^a m - 

a (z, z)O m (0) ~ 4 (3C W - Q - 40 A ) 

zz 

The resulting /5eta functions are 

/9 a = -2fi£ + 2# a # m + g A (g a + g m ) 

(3 m = -2g^ - Qg a g m + 3g A (g a + g m ) (3.10) 
Pa = 2g A + Sg a (g a + g m ) 

If one starts with a model with only the random spin potential a, i.e. only g a ^ 0, 
then g A and # m are generated under renormalization. Note that the /3eta functions simplify 
dramatically when g a + g m = 0, which is the subject of the next section. We will return 
to a general analysis of the possible universality classes contained in these /3eta functions in 
section V. 
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IV. SPIN-CHARGE SEPARATION 



A. Symmetries of the effective action 

The free conformal action S c f t has a maximal osp(4|4)i current algebra symmetry. The 
conformal currents J A generating this symmetry correspond to all possible bilinears in 
fermions and ghosts. Suppressing the possible index structure, J A = {ipip,ipf3, (3/3}. The 
perturbing operators are bilinear in these currents: O a = d^B J A J for some tensor d^B 
and similarly for O m and Oa- 

For arbitrary g a ,g m ,gA, the osp(A\4) symmetry of S c f t is broken but an osp{2\2) sym- 
metry is preserved. To describe the symmetry, let us define the left-moving currents: 



h = y,f+P-> s ± = $>i/?; (4.i) 

i i 

where = —eji, €12 = 1, and similarly for the right-movers J = ij) + ip_,H = /3,/3_, 

These currents generate an osp(2\2)k current algebra at level k = —2. The non-zero operator 
products are 

J(z) J(0) ~ -- 2 , H(z)H(0) ~ - 2 

Z A 

2 2k A 

J(z)J±(0) ~ ±- J±, J + (*)J_(0) ~ — J 

z z z z 

J(z)S±{Q) ~ ±-S±, J{z)S±(0) ~ ±-S± 

z z 

H(z)S±(0) ~ ±^S±, H{z)S±(0) ~ ^-S± 

J±(z)S T (0) ~ -£±, J±(^)5 T (0) ~ --S± 

£±(^£±(0) ~ ±- z J± (4.2) 

5 + (^)5_(0)~4 + i(iJ- J) 

5 + (z)iS_(0) ~ - + - (if + J) 

The currents J and J± generate a charged su(2) subalgebra under which the fermions if)± 
transform as doublets. 

The osp(2\2) symmetry is present for arbitrary potentials before disorder averaging, as 
we now describe. The fermionic generators define nilpotent transformations of the bosonic 
and fermionic fields. The action ( |2.9|) is then an exact variation with respect to these 
transformations, making then clear the fact that this osp(2\2) symmetry is preserved. For 
example the generator S+ induces the transformation 5ip l + = 0, 8ip % _ = PL, S/3 l + = 
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8f3 l _ = 0. The symmetry acts left-right diagonally, so for the right movers 8ip + = 0, 8ip_ 
~/3 l _, 8(3, = 8~0 l _ = 0. The action (j2.9|) may be written as 



where 

9 = - !p_p + ) + id ■ (ip-a]3 + + ~$_o(3 + ) + %A ■ (tp-a(3 + + ij>_aj3 + ) (4.3) 

Hence 8S = since 8S c ft = and 8 2 = 0. This holds similarly, but with another operator 
©, for £_ which generates the transformations 8ip l + = e^(3 J _, 8ip l _ = and 8f3 l + = e^ipL, 
8f3 l _ = 0. This osp(2\2) symmetry would not be preserved if we had chosen the random 
gauge potential in u(2) or if had added extra scalar potential randomness. 

This implies that after disorder averaging the perturbations by O a , O m and Oa in 
the effective theory preserve a global, left-right diagonal, osp(2\2) symmetry. The current 
conservation law takes the left-right diagonal form 

djJ a + d z T = (4.4) 

for J a , J a any of the eight osp(2|2)_ 2 currents. 

Let us describe more explicitly how the osp(2\2) symmetry is manifested in the effective 
theory. In the sequel we will set g a = —g m = 9 and the effective action will contain the 
operator O g = O a — O m . Many terms cancel in the combination O a — O m and the result 
can be written as an osp(2\2) current-current perturbation. This simplification is analogous 
to the gv + Qm = line for the random U(l) fermions which has gl(l\l) symmetry, studied 
T0|| • By repeated use of the identity 



m 



a ij a nm = 2$im8jn ~ 8ij8 nm (4.5) 

one can express O g in terms of the above osp(2\2) currents: 

O g = O a - O m (4.6) 
= - J J + HE + hj-J+ + J+1-) + SJS + - S + S_ - S_~S+ + S+S- 

The operator O g has the structure of the quadratic Casimir for osp(2\2). Namely, O g = 

J2 a ,bCabJ a J b , where J a are osp{2\2) currents and C a b corresponds to the quadratic Casimir. 
Thus O g is osp(2\2) invariant. 

The operator O a + O m cannot be written only in terms of the osp{2\2) currents; one 
needs some of the osp(4|4) currents. Let 
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u% = W±P± + Al> i ±) ; ?- = iMftf ; v+ = P-*i>+ (4.7) 



Then, 



O a + O m = -Og + O 



s 



6 = \ (r/'ir-! + /r/77") + ~ (t f ■ L f - i h ■ u 



with 



+U 1 JU+ - U%U_ + - \V+ ■ V- - V- ■ V+J (4.8) 

We have already seen that O g is osp(2\2) invariant. One can check explicitly that O is also 
osp{2\2) invariant. 

Let us turn now to the perturbation O 'a- This operator as defined is a left-right current- 
current perturbation: 

A = j2L a T (4.9) 

a=l 

where 

L a = L a f + L a b = V>_ct> + + [3-a a f3 + (4.10) 

The central extension (level) cancels between the fermions and ghosts and the result is that 
L a generate an su{2)^ current algebra at level k = 0: 

L a (z)L b (0) ~ -J abc L c {Q) + reg. (4.11) 

Again, since Oa takes the form of the Casimir of su(2), Oa preserves a global su{2) sym- 
metry. Furthermore it is easy to check that this global su{2) commutes with the osp(2\2) 
symmetry: 

L°, J b ] = 

where J a are the osp(2\2) currents defined in eq. ( |4.1| ). 

In summary, the model has a global osp{2\2) ® su{2) symmetry for general g a ,g m ,gA- 



B. Quasi-spin-charge separation in the conformal field theory 

Since O g only involves the osp(2|2)_ 2 currents and Oa the su(2) currents, it is important 
to understand how the conformal field theory decomposes in terms of these current algebras. 
It is straightforward to check that the su(2) and osp{2\2)_ 2 current algebras commute, so 
that the conformal field theory defined by S c f t in eq. ( |3.3| ) has an osp(2|2)_ 2 <8> su(2) 
symmetry. 

We now show the much stronger result that the full stress tensor, which is the Sugawara 
stress tensor for the osp(4|4)i supercurrent algebra, separates into two commuting pieces: 

T os p{A\i)i = -^osp(2|2)_ 2 + T su (2) (4-12) 

All T's in the above equation are Sugawara stress tensors and have Virasoro central charge 
equal to zero. The above equation is proven in the appendix. 

Since osp{2\2)^ 2 contains charged currents whereas su(2) does not, the equation (|4.12| ) 
implies a kind of spin-charge separation and we will use this terminology in the sequel. 
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Based on the separation ( |4.12| ) one would expect that the Hilbert space of osp(4|4)i 



factorizes. However by trying to perform explicitly this factorization in some simple cases, 
we found that it is not possible: 

^osp(4|4)i 7^ ^osp(2|2)_ 2 ® 7~tsu{2) 

This non-factorization is intimitely related to the way the fields that appear in O transform 
under osp(2\2). This is described in Figure 2. In particular, the su(2) currents L a turn out 
to be susy exact, e.g.: 

S + (z)V*(0) ~ -L a (4.13) 

z 

Let us argue for the non-factorazibility of the Hilbert space by contradiction. The 24 fields 
of O form three copies of an eight dimensional reducible but indecomposable representation 
of osp(2\2), that we shall denote [8]. Consider the state \L a A ) = {L% x — L^-JIO), where 
Lti = § dz L a (z)/2mz. If the factorization of "Hosp^^-i holds, then 

\L\) = E \° s P)i ® \su)i G [8] ® [3] 

i 

where [3] is the adjoint representation of su{2). Let us now act on |L^) with the su(2) 
current. As a state in 7^ sp(4|4)i one finds that L^ 1 |L^) = 5 ab |0F OC k) with |0F OC k) the vacuum 
of the Fock space. As a state in r H osv {2\i)^ 1 ® 7^su(2) > one would find 

£-il^> = E \° s v)i ® L\\su)i E [8] ® [1] 
where [1] is the su(2) singlet. Thus if the factorization holds, one would deduce that 

|0 F ock) e [8] ® [l] 

Since |0pock) is annihilated by all osp(2|2) generators, this would mean that |0p O ck) is in 
the image of susy generators, i.e. there would exist a state \Q) such that for example 
|0Fock) = iS-lfi). This is clearly a contradiction. A similar argument shows that eq. ([4.13 ) 
leads to a contradiction in the factorization of V^. 

Another argument for the non-factorization is based on the fact that the current alge- 



bras osp(2|2)_2 and su(2)o lead to logarithmic correlation functions (13| ||14|| , whereas the 
osp(A\A)i is a free theory with no logarithms. It appears impossible for these logarithms to 
cancel in the product of two logarithmic correlation functions. The logarithmic nature of 



these theories can in fact be traced to transformations such as ( 4.13Q fT3f 



We now describe a few of the results we need concerning these current algebra conformal 
field theories. Representations of osp(2\2) are characterized by the quantum numbers of the 
su{2) ® u(l) subalgebra generated by (J, J±) and H. Highest weights can be labeled (j, b) 
where j = 0, 1/2, 1, .. is the spin of the charge-su(2) and b = H/2. These representations have 
dimension 87. The conformal scaling dimension (left or right-moving) of the corresponding 
primary fields, as determined from the Sugarawa construction, is 
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On the other hand, the primary fields of the su(2)k current algebra are characterized by the 



spin j of the spin-sw(2) only and have conformal dimension [] I 



A- 



«*(2) _ j{j + 1) 



k + 2 



(4.15) 



A simple check of eq. ( |4.12| ) is based on the dimensions of the eight original fermion and 
ghost fields ip±, fi±- in the original conformal field theory these fields have conformal dimen- 
sion A = 1/2. Under the osp(2\2) they transform according to the two four- dimensional 
representations fit, and (^i, ip+, P+, P-) corresponding to (j = 1/2,6 = 0) and 

its conjugate. At level k = —2, these have dimension A osp ( 2 ' 2 ) = 1/8. The same fields trans- 
form as spin j = 1/2 doublets according to the su(2) , and have dimension A su ^ = 3/8 
at level zero. The exact decomposition ( 4.12| ) implies that these dimensions must add up 
properly: 1/2 = 1/8 + 3/8. 



C. Infrared fixed points and critical exponents 

We now study the model when g a + g m = 0. Setting g a = —g m = g, the effective 
action contains the current-current operators O g and O a . The operator O which couples 
the osp(2|2)_ 2 an d su(2) current algebras is not present. The model in || was mapped 
onto a super-spin chain with osp(2\2) symmetry and Heisenberg-type hamiltonian. This 
further supports the identification of this super spin chain with our model on the line g a + 
g m = 0, since, at least for ordinary bosonic algebras, Heisenberg hamiltonians correspond 
to symmetry preserving current-current perturbations built on the quadratic casimir in the 
continuum limit. (As Lie super algebras, osp(2\2) and s/(2|l) are identical.) 

The effective action now contains interactions that do not couple the spin and charge 
currents: 

S eff = S cit + ( (gO g + g A A ) (4.16) 

J Z7T 

This decoupling is of course consistent with the more general /Seta functions (|3.10|) . Setting 
g a = -g m = 9 one finds 

(3 9 = -Ag\ (3 gA =2g\ (4.17) 

As explained in section II, the model of Senthil et. al. corresponds to g and g A positive. 
From the above /3eta functions one sees that g A is marginally relevant whereas g is marginally 
irrelevant. In order to deduce something exact concerning the infrared (IR) fixed point theory 
from the one-loop /3eta functions, one needs to make a hypothesis concerning the role of the 
higher loop corrections. Consider for comparison the su(N) Gross-Neveu models, which are 
su(N) current-current perturbations, with /3eta functions as in eq. ( |4.17|) , i.e. (3 g = g 2 in 



a certain convention. When g > the perturbation is marginally relevant, i.e. g grows at 
large distances. Higher loop corrections do not modify this, i.e. the flow to the IR does not 
stop at some finite value of g corresponding to a non-trivial fixed point. Rather, g eventually 
flows to infinity. The theory is thus massive, and in the infrared all these massive modes 
disappear leaving no massless degrees of freedom, i.e. the infrared theory is an empty theory. 
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When g < 0, the perturbation is marginally irrelevant, i.e. g flows back to zero and the 
unperturbed conformal current algebra is recovered in the IR. We will make the hypothesis 
that this is the only possible behavior for general current-current perturbations, i.e. the 
only possible fixed points are g = or oo. 

For our model we then have the following picture. In the IR, gA flows to infinity and 
the spin degrees of freedom are massive and decouple. The coupling g on the other hand 
is marginally irrelevant. Therefore the IR fixed point is the coset osp(4|4)i/su(2)o, and the 
theory arrives in the IR via the operator O g . Due to the non-factorizability of the Hilbert 
space, this coset conformal field theory is not precisely the osp(2\2)_ 2 current algebra, even 
though it possesses this current algebra as a symmetry and the conformal dimensions of the 
coset theory are the same as the osp{2\2) current algebra because of (|4.12j ). 



This identification of the IR fixed point allows the computation of certain critical expo- 
nents. The density of states p(E) is 

p(E) = -^Tr 8{H — E) = \ lim lmTY 1 . (4.18) 

V 7lV e^0+ H — E — IE 

where V is the two-dimensional volume. This implies that the disorder averaged density of 
states is proportional to the one-point correlation function: 



p(E) oc Re($ E ) (4.19) 



Let equal the scaling dimension of $g. Then, since the action (|2.9|) is dimensionless, 
viewing £ as a coupling, dim(E) = 2 — T E . Since E is the only dimensionful coupling in 
the theory one deduces 

pjE) (x E Te/{2 ' Ve) as E^O (4.20) 

One can also define a correlation length 

CeocE-"*, v E = l/{2-T E ) (4.21) 

(In ||, u E was referred to as z/g.) 

For our theory, = + ... and T E is the scaling dimension in the IR, which follows 

from the osp(2|2)_ 2 conformal dimension of ijj± } which as explained above is A/^o? = V^- 
Thus T E = 1/4, and: 

pjE) oc E 1/7 , v E = 4/7 (4.22) 

Numerical simulations of the osp{2\2) invariant spin-chain agree very well with v E = 4/7 
0. Since the system flows toward the infrared fixed point along a marginal direction, the 
scaling (|4.22|) is up to computable logarithmic corrections. 

In the above scenario, the su(2) spin degrees of freedom are becoming localized, hence 
the terminology "spin quantum hall effect" . The above model can also exhibit charge local- 
ization by changing the signs of the couplings. Namely, if g and gA are negative, then g is 
marginally relevant, i.e. it flows to — oo in the IR and gA is marginally irrelevant. Following 
the same reasoning as above, in this case the charge degrees of freedom decouple in the flow. 
The IR fixed point is now the coset osp(4|4)i/ osp(2|2) _2- Here, T E = 3/4, and: 
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p(E) oc E 3 / 5 , 



£e oc E ' for charge delocalization 



(4.23) 



Finally, if g < 0, > 0, then g, g^ flow to — oo and oo respectively, and both sectors are 
massive and decoupled in the IR. This implies = 0, and a constant density of states at 
E = 0. 

D. Path Integral Factorization 

Let us now present a path integral derivation of the spin-charge separation which has 
the advantage of being non perturbative and thus valid all along the RG trajectory. One 
can view this formulation as a way of defining the coset osp(4\4)i/ su(2)q. It consists in 
decoupling the random gauge field A using chiral gauge transformations. The spin-charge 
separation will then be a consequence of the fact that the osp(2\2) currents are invariant 
under the chiral gauge transformations. This decoupling is similar to the solution of the 
random gauge potential described in ||. 

Let S be the fermionic action ( |2.9| ) and Z(A,a,m) its partition function Z(A,a,m) = 
J D*5l e~ s ^\ At a fixed realization of disorder, the gauge potential A may be gauged away 
by a chiral gauge transformation by parameterizing A as: 

iA = G~ 1 d s G ; iA = G*d z G*- 1 (4.24) 

with G an element of the complex su(2) c group, i.e. G is a two by two complex matrix with 
determinant one. This is always possible on the sphere. This parametization is such that 
+ iA)ip + = (ip^G~ 1 )dz{Gip+)- Let us now denote by ip' the chiral gauge-transformed 
fermions: 

= ip_G~ 1 ; = ^_G* (4.25) 
1/+ = Gil> + ; % = G*~ X ^ + 

Gauge transformed bosons f3' are defined similarly. The fermionic action ( |2.9| ) may be 
rewritten in terms of G and . It becomes: 

72 

+ i^G(a-a-m)G*^' + + 0_G*- 1 {a-a + m)G- 1 'ilj' + ) (4.26) 

This is the effective action for the Dirac operator coupled to the disorder variables Ga ■ aG*, 
G* ~ 1 a-aG~ 1 and GmG*. As we shall see below, the field G will be coupled to the spin degrees 
of freedom whereas the fermions shall be coupled to the charged sector. The density of 
states has a factorized expression in terms of these new spin and charged variables: 

p{E) oc Re (GG*) (4.27) 

We must not neglect to take into account the jacobians of the transformations A — > G 
and \& — > as well as to represent the partition function Z(A, a, m). The jacobians, which 
are determinants of Dirac operators, are computed using the chiral anomaly. They may be 
expressed with the help of the WZW action: 
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DA 



DG 



Det (i0 + A)= exp (S WZW (GG*) ) 
exp (AS WZW (GG*) ) 



with <9 tUifu , the WZW action. 

To represent the partition function, note that Det (i0 + A) = Z(A,a = 0,m 
Therefore, chiral gauge transformations applied to the bosonic beta system give: 



0). 



Det (i0 + A) 



Z(A, a, m) 

Gathering the jacobians we get the action 



J Dl3'exp( -T(/3'\G,a,m) 



S = -AS WZW (GG*) + T(¥\G, a, m) + T{fl\G, 



(4.28) 



(4.29) 



with T defined in eq.(4.26). This is the action to compute correlations using G, ty' and (5' 
as the path integral variables at fixed disorder. One still has to add the disorder measure 
( [2~3l , 12.5|) to compute averaged correlations. 

In eq.( [4.29| ) the random spin variables A or G are not yet decoupled. This decoupling 
only appears when g a + g m = 0. Indeed, in that case integrating over the disorder a and m 
yields to current- current type interactions as in eq. Q4.6D but with osp(2\2) currents bilinear 



in the fermions $ instead of the fermions iff'. The crucial point is now to remark that these 
osp(2\2) currents are invariant under the chiral gauge transformations ( |4.25|) . For example, 



GG' 1 ip. 



m G inG jm ^ij 



J- 



where we used that G in Gj m 6ij 



detG e nm = e nm since G has unit determinant. Similarly, it 
is easily checked that all osp(2\2) currents are invariant under chiral gauge transformations. 
Note that this is true because the spin disorder variables belong to su{2). This would not 
be valid if for example the spin disorder variables were taking values in u(2) instead. 

Hence, after integrating over the disorder at g a + g m = 0, the spin random variables G 
decouple from the and f3' system. This is the spin-charge separation. We are thus left 
with the action: 



S = -AS WZW (GG*) + r(tf'|G = 1, a, m) + T((3'\G = 1, a, m) 



(4.30) 

with the same disorder measure for G and a, m as in eq. ( pT3|j2~5| ) with g& arbitrary but 
g a + g m = 0. The first term describes a WZW theory on the coset space su{2) c / su(2) at 
level k = —4. This may be thought of as the theory 'inverse' to the su{2) WZW theory at 
level k = [ 15| , since the conformal dimensions of primary field in the two theories have 
opposite sign. The second and third terms describe the osp(2\2) current- current perturbation 
of the ty' — f3' system. As in previous section this ^' — (3' system may also be described as 
an osp(A\A) WZW model at level one and the current-current perturbation only couples to 
the sub-sector generated by the osp(2\2) currents. 

In the last section we argued that g^ flows to infinity under RG flow. At g& = oo, the 
measure ( [2.3| ) on G is flat and only the WZW action S WZW (GG*) at level k = —4 remains. 
This describes the spin sector. 
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If g a = —g m > 0, the osp(2\2) current perturbation is irrelevant and we recover the 
previous description. Recall the original fermion from which the density of states is 
computed, are related to G and ^' by the chiral gauge transformation ( [4.25| ), so that the 
density of states is factorized as in eq.( |4.27| ). 

If g a = —g m < 0, the current- current perturbation is relevant. We may then propose that 
in the infrared the ^' — (3' system is described by the coset theory osp(A\A)k=i/ 'osp(2\2)k=-2, 
which is equivalent to an su(2) model at level zero. As a consequence, the fermions shall 
flow in the infrared to fields with scaling dimension opposite to that of G, and the original 
fermions \1/ shall flow to fields with zero scaling dimension. This suggests that in this case 
the density of states is finite and regular at zero energy. 



V. RG PHASE DIAGRAM 

In this section, we return to the general model with g a + g m 7^ 0, and describe the global 



features of the phase diagram based on the one- loop RG equations (|3.10|) . 

Our method consists in extracting the asymptotics of the RG trajectories by looking for 
directions in the coupling constant space which are preserved by the RG flow. Then, to 
analyze whether these asymptotic trajectories are attractive or not, stable or unstable, we 
project the RG flow onto the sphere and study the vector field thus obtained. This will allow 
us to point out the special role played by perturbations along a so-called strange direction. 

Let us first look for lines in the coupling constant space, with coordinates g = (gA, g a , g m ), 
which are preserved by the RG flow and which pass through the origin. These correspond to 
trajectories which are straight lines and therefore for which the RG velocity field g = (3(g) 
is co-linear to the vector g. The equations for these fixed line trajectories are thus: 

P(g)Ag = (5.1) 

The above equation is equivalent to finding solutions of the RG equations of the form 
g l (t) = x l X(t) where x % are constants independent of the RG time t = log I and d\(t) jdt = A 2 . 



Substituting this into (|3]6|), one finds 



x — CjfcX^x (5.2) 

where as before Cj fc are the operator product coefficients. With the linear relation g % = x l \ 
among the couplings, the effective action contains a single running coupling constant: 

72 

S cS = S cft + A J O x , O x = J2 x'Oi = * A Oa + x a O a + x m O m (5.3) 



For x % a solution to ( [5. 2D , the operator O x closes on itself under operator product expansion: 

O x (z,z)O x (0) ~ 0,(0) (5.4) 

To lowest order the /3eta function for this specific perturbation is (3\ = A 2 + • • • 
There are six solutions to eqs.( |5.1|) (or eqs. (|5.2|) ). They are: 
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A= (g A ^ 0, g a = 0, g m = 0) (5.5) 

B= (g A = 0, g a = 0, g m + 0) (5.6) 

C = {g A = 0, g a + g m = 0) (5.7) 

D= (g a + g m = 0, g A + 2g a = 0) (5.8) 

E= (6^ a = -3(l+V2)^A = -2(l + 2v^)p m ) (5.9) 

F= (6^ = 3(^-1)^ = 2(2^-1)0 (5.10) 



The two first lines correspond to the simple models with only random gauge potential |T6|,|6[| 



or only random mass |17Hl9i. The third and fourth ones correspond to the case we discuss 
at length in section IV. Contrary to the first four solutions, the last two do not seem to have 
an obvious simple algebraic interpretation 0. 

These six solutions correspond to twelve possible stable directions because we have to 
choose an orientation on the line. Only three of them are in the domain of positive coupling 
constants g A > 0, g a > and g m > 0. The first two, which are the solution A with g A > 
and the solution B with g m > are in the border of that domain. The third one, which is 
the solution F with g a > 0, sits in the middle of the domain of positive g's. We shall call 
this solution the strange direction. 

To analyze deeper the RG flow let us now project it onto the sphere. This is worth 
doing since the beta functions are homogeneous. Since there are three coupling constants, 
we may parameterize them with two angles, that we shall denote 9\ and 62, and the radial 
coordinates p, p 2 = g\ + g 2 ^ + g^. The RG equations ( p.!0|) may then be written as two 



equations for the angular variables 

^ = p/3;(0l,0 a ) (5-11) 
together with one equation for the radial variable: 

p = p 2 2 ) 

The explicit expressions for the vector fields (3j or (3 P are easy to find. In the angular 



equations (|5.11 ) we may absorb the factor p into a redefinition of the parametrization of 
the RG trajectories. This does not change the topology of the RG curves but only the 
speed at which the RG trajectories flow on these curves. So we shall analyze the vector field 
9j = i3j(9i,9 2 ) on the sphere. See figure 1. 

The vector field j3j has twelve zeroes which correspond to the twelve fixed directions 
( |5.5| - |5.10| ). We can compute the sign of p in these directions to know whether these straight 
RG trajectories are escaping to infinity (p > 0) or are flowing back to the origin (p < 0). 
To decipher whether the fixed directions are attractive or not we have to analyze whether 
the corresponding zero of the vector field f3j on the sphere is attractive or not. For that we 
linearize the vector field (3j on the sphere at its zeroes and compute its eigenvalues. Positive 



2 The operator O x is bilinear in the currents, O x = J A d A sJ with J A the osp(4|4) currents and 
d A B some bilinear form. The condition ( |5.4j) is then: ff B d A jqdBMfK M = diK- The algebraic 
interpretation of that equation is not clear to us. 
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eigenvalues correspond to repulsive fixed directions, zero eigenvalues to locally marginal 
directions. For a fixed direction to be generically attractive, the two eigenvalues have to be 
non positive; otherwise to be attracted to the fixed direction requires fine tuning. The result 
is: 
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The indices ± refer to trajectories flowing to infinity or back to the origin. 

There are only four fixed directions which are generically attractive: they stand on 
the directions B + , C+, D- and on the strange direction F + . They all are asymptotes to 
RG trajectories flowing to infinity except LL which corresponds to asymptotic direction of 
trajectories looping back to the origin. 

Let us pause to reconsider the case g a +g m = is this language. In that case the solutions 
of the RG equations are 1/gA — ^-/9a — and l/g a — 1/ — At. If g A > and > 0, 
then gA increases whereas g a decreases. These trajectories are asymptotic to the direction 
A since gA blows up at a finite value of t at which g a is finite. However the blow up times 
are not physical since for qa large enough the one-loop analysis is no more valid and one has 
to relies on the non-perturbative analysis done in previous section. If g\ > but < 0, 
then both couplings increase in magnitude. Which one blows up first depends on the initial 
data. The asymptotic trajectories are then either the direction A or C depending if g A + 2g^, 
is positive or negative, so that the line gA + 2g Q = is a separatrix in the phase diagram. 
On contrary, ii g A < but g® > 0, the couplings decrease and flow back towards the origin 
along the direction D. 

The strange direction F + is the asymptotic direction for all trajectories starting initially 
with positive coupling constants g a , m ,A, except those which are fine tuned to be in the 
direction A or B. The domain gj > is stable under one-loop RG: no trajectories can 
escape from it. The fact that there is one and only one asymptotic direction in the domain 
of positive coupling constants reflects the universality of that behavior: whatever the values 
of the initial coupling constants, the system flows along this direction in the infrared. The 
strange direction is hence the one which should be relevant to the description of the low 
energy behavior of the su(2) random Dirac operators, e.g. the network model described in 
section II. Furthermore, since p > for the strange direction, it is a strongly coupled system. 
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In summary, based on the one-loop RG equations, it appears the generic network model 
is in a different universality class than the one fine-tuned to the line g a + g m — 0. However it 
certainly remains a possibility that the higher loop corrections can spoil the above analysis. 



VI. DISCUSSION 

In summary, we have shown that the network model with the identification of couplings 
9a + 9m = exhibits a spin-charge separation in the effective disorder averaged theory, and 
this allows a precise identification of the infrared fixed point as the coset osp(4|4)x/ su(2) . 
This coset conformal field theory has some novel features in that it possesses an osp(2|2)_2 
supercurrent algebra symmetry, but due to the non-factorizability of the Hilbert space, it 
is not identical to the osp{2\2) theory defined by the Sugawara construction. The resulting 
critical exponents agree with the predictions based on percolation j|] and the numerical 
simulations of the super spin chain studied in M. Based on the one-loop renormalization 
group, we have argued that the network model without the constraint g a + g m = is in a 
different universality class. 

How our analysis relates to percolation remains an interesting open question. There 
appear to be two possibilities. One is that the 1-loop strange direction we described in 
section V survives to higher loops and indeed corresponds to percolation; our coset the- 
ory osp(A\A)i/ su(2)o is then the wrong fixed point. Another possibility is that the higher 
loop corrections actually restore the symmetry to the line g a + g m = 0, and the coset 
osp(4\4)i/ su(2)o is a new description of percolation with osp(2\2) symmetry. The latter 
possibility can be investigated by comparing the conformal blocks of percolation with those 



computed in [20 



It would be interesting to construct explicitly the continuum field theory corresponding 
to the super spin chain along the lines of [^TJ. It is important to understand the role of the 



global osp(2\2) symmetry of the spin chain if this continuum limit indeed corresponds to 
percolation. 
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VIII. APPENDIX: SPIN-CHARGE SEPARATION OF THE STRESS TENSOR. 



We give here some details on the proof of eq. ([4.12|) for the energy momentum tensor. 



Recall that the Sugawara construction of the energy momentum tensor in terms of the 
currents J a is: 

/ fi Sjab \ 

T(z) — k lim I J a (w)C ab J b (z) - k (8.1) 

w^z \ ( z - w )2 
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with C a b the quadratic Casimir and the normalization constant k is chosen such that the 
currents have conformal dimension one. In practice T is computed by extracting the regular 
term in product J a (z)C ao J b (0). 

For the osp(2\2) algebra at level —2 this gives: 

T osp(2 | 2) „ 2 = l -(j 2 -H 2 - \{J-J+ + J+J-) + (S + S_ - S_S + ) + (S_S + - S + S_) ) 

The osp(2\2) currents are given in eq.( [4.1|) in terms of the (3 — if) system. This Sugawara 
energy momentum tensor can thus be expressed in terms of these bosonic and fermionic 
fields. A simple computation yields: 

T„ SP ( 2 |2)„ 2 = \{d z ifj.ifj + + d z ip + if)_ + d z (3.f3 + - d z [3 + P_ - 4ftM + )(/M+) 
+ 3(^ + ) 2 - (P.p + ) 2 + 2(^+)0M+)) 



Similarly, with the normalization (57?) for the Pauli matrices, the su(2) Sugawara tensor 

is: 

Ti r a r a 

su(2) — g L> Li 

with L a defined in eq.( [4.10[ ). Again with the help of the identity ( |4.5| ), this may be written 
in terms of the (3 — if) system as: 



T, 



su(2) = ^(m^+ + d z i/>+i/>- + d z p_P + - d z p + P_) + A(ip_f3 + )(f3_ifj + ) 
- 3(^+) 2 + {(3.(3 + ) 2 - 2ftM + )(/3-/3+)) 
Adding both pieces we get: 

r osp (2|2)_ 2 + T su ( 2 )o = g ( d zi>-^+ + d z ip+if)- + d z (3-(3+ - d z (3 + f3-) 
This is the energy momentum tensor T for the conformal field theory S c ft- 
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FIGURES 

FIG. 1. The vector field [3j on the sphere. We choose the coordinates ^g a = cost/, 
2(g a + g m ) = sinqcosp and 2gA = singsinp with q 6 [0, 7r] and p G [0, 27r]. 
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FIG. 2. Structure of the indecomposable osp(2\2) representation formed by the field appearing 
in O, Ls = Lj + Li, = L and La = Lj — L^. The solid arrows describe the action of S- and the 
dashed ones the action of S- . 
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